Mon. Not. R. Astron. Soc. 000, 000-000 (0000) Printed 19 July 2011 (MN I*TeX style file v2.2) 

The relation between metallicity, stellar mass and star 
formation in galaxies: an analysis of observational and 
model data 



o 



Robert M. Yates^*, Guinevere Kauffmann^ & Qi Guo^ 

^ Max Planck Institut fur Astrophysik, Karl-Schwarzschild-Str. 1, 85741, Garching, Germany 
National Astronomical Observatories, Chinese Academy of Sciences, Beijing, 100012, China 



vn 



Accepted ??. Received ??; in original form 



o 
u 

6 



> 

in 

m 

o 



% 



ABSTRACT 

We study relations between stellar mass, star formation and gas-phase metallicity in 
a sample of 177,071 unique emission line galaxies from the SDSS-DR7, as well as 
in a sample of 43,767 star forming galaxies at 2: = from the cosmological semi- 
analytic model L- Galaxies. We demonstrate that metallicity is dependent on star 
formation rate at fixed mass, but that the trend is opposite for low and for high stellar 
mass galaxies. Low-mass galaxies that are actively forming stars are more metal-poor 
than quiescent low-mass galaxies. High-mass galaxies, on the other hand, have lower 
gas-phase metallicities if their star formation rates are small. Remarkably, the same 
trends are found for our sample of model galaxies. By examining the evolution of the 
stellar component, gas and metals as a function of time in these galaxies, we gain 
some insight into the physical processes that may be responsible for these trends. We 
find that massive galaxies with low gas-phase metallicities have undergone a gas-rich 
merger in the past, inducing a starburst which exhausted their cold gas reservoirs and 
shut down star formation. Thereafter, these galaxies were able to accrete metal-poor 
gas, but this gas remained at too low a density to form stars efficiently. This led to a 
gradual dilution in the gas-phase metallicities of these systems over time. These model 
galaxies are predicted to have lower-than-average gas-to-stellar mass ratios and higher- 
than-average central black hole masses. We use our observational sample to confirm 
that real massive galaxies with low gas-phase metallicities also have very massive black 
holes. We propose that accretion may therefore play a significant role in regulating 
the gas-phase metallicities of present-day massive galaxies. 

Key words: ISM: abundances - Galaxies: abundances - Galaxies: star formation - 
Galaxies: fundamental parameters - Galaxies: evolution 



1 INTRODUCTION 

iVietals are ubiquitous througliout galaxies. They are syn- 
thesised in stars and liberated into the interstellar medium 
(ISM) when stars shed their outer gaseous envelopes towards 
the end of their lives, and in some cases also into the in- 
tergalactic medium (IGM) when the highest mass stars ex- 
plode as supernovae. The amount of metals in the diffuse 
gas around galaxies also determines the rate at which it is 
able to cool and form stars. Metallicity is, therefore, one of 
the key physical properties of galaxies, and understanding 
the processes that regulate the exchange of metals between 
stars, cold interstellar gas, and diffuse surrounding gas can 
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help us understand the physical processes that govern galaxy 
evolution in general. 

The metallicity of stars and gas in galaxies is known 
to correlate strongly with the ir luminosities, circu l ar veloc- 
ities and stellar masses (e. g. Leaueux et all 19791 : iGarnetj 
l2002l : iTremonti erall l2004l : TGallazzi et al.ll2005l '). However, 
the physical processes that drive these correlations are not 
yet fully understood. 

iMathews fc Baked (|l97ll ) and Ibarsonl (|l974l ') first sug- 
gested that interstellar gas can be driven out of galaxies by 
supernova explosions as galactic outflows. They predicted 
that galaxies of smaller mass have lower metal abundances 
because their lower escape velocities allow freshly enriched 
gas to be more efficiently removed. 

The input of energy from supernova explosions is now 
routinely incorporated into hydrodynamical simulations of 
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galaxy formation, either in the form of thermal heating or 
'kinetic feedback', whereby radial momentum kicks are im- 
parted to particles surrounding sites of star formation in the 
galaxy. Although these simulations are now able to demon- 
strate that galactic outflows can yield a good match to 
the observed mass-metallicity relation (e.g. iKobavashi et al.l 
120071 : IScannapieco et aLlbOOa ). this does not mean that out- 
flows are the only process at work in regulating the metal- 
licit ies of real galaxies. 

iBrooks et al.l ((20071) argue that a mass- dependent star 
formation efficiency (SFE) is required to explain observa- 
tions, in addition to supernova feedback. In this scenario, 
less massive galaxies convert their gas reservoirs into stars 
over longer timescales than more massive galaxies. There- 
fore, less massive galaxies have higher gas-to-stellar mass 
ratios and are consequently less metal-rich. Motivated by 
their work on momentum-driven wi nds in smoothed particle 
hydrodynamics (SPH) simulations, iFinlator fc Davg (|2008l ) 
have also suggested that both metal-rich outflows at all 
masses and a variable star formation efficiency must play 
roles in explaining the observed mass-metallicity relation for 
nearby galaxies. 

Alternatively, iDalcanton et al.l (1200J) suggest that 
metal-poor infall can regulate metallicity in disc galaxies. 
Given that lower mass galaxies tend to have lower star for- 
mation rates (due to the fact that their cold gas densities 
are lower), a net dilution takes place when the time-scale 
for star formation falls below the time-scale for accretion of 
metal-poor gas. This would drive down the low mass end of 
the mass-metallicity relation. 

Finally, variations in the initial mass function (IMF) 
have also been cited as anothe r factor that might influence 
the m ass-metallicity relation. iKoppen. Weidner fc Kroupal 
112003) propose that a SFR-dependent (and therefore stellar 
mass-dependent) IMF causes different galaxies to produce 
different effective oxygen yields. This hypothesis is based 
on the premise that most stars form in stellar clusters and 
that smaller, less actively star-forming galaxies are domi- 
nated by clusters with lower masses, containing a smaller 
fraction of massive, oxygen-producing stars. We note, how- 
ever, that there is little observational evidence for system- 
atic IMF variations between stellar clusters ( e.g. lElmegreenI 
I2OO6I : iFumagalli. Da Silva fc KrumholjIJoTil ') . 

The fact that the interpretation of the mass-metallicity 
relation is subject to considerable ambiguity has prompted a 
number of authors to consider alternative ways of quantify- 
ing metallicity in galaxies. For example, higher dimensional 
relations that include additional physical properties could 
provide better const raints on the processe s that regulate 
metallicity. Recently. iMannucci et al.l (|2010f ) have proposed 
such an extension to the mass-metallicity relation, called the 
fundamental metallicity relation (FMR). In that work, it was 
shown that the gas-phase metallicities [Z) of both local and 
high-redshift galaxies were dependent on both stellar mass 
(Mt) and star formation rate (SFR). The FMR provides a 
prediction of the metallicity of local galaxies with a fa scat- 
ter of only ~ 0.05 dex. This is a substantial improvement on 
the m ean scatter of ~ 0.1 dex reported by iTrcrnonti ct al. 
l|2004l for the M,-Z relation. Metallicity was also found to 
be strongly dependent on SFR at low stellar masses, but 
onl y very weakly dependen t on SFR at high stellar masses. 

iMannucci et all (|201(]| ') further showed that the FMR 



describes galaxies out to z ~ 2.5. The observed evolution of 
the Mt-Z relation was therefore attributed to migration of 
galaxies along the FMR plane to higher masses and lower 
star formation rates over time. 

Prompted by these findings, in this work we study 
higher dimensional relations between metallicity and a va- 
riety of physical galactic properties. We examine whether 
the Mt-Z relation exhibits additional dependences on SFR, 
specific SFR (sSFR) and Afgas/Af, . Comparisons are made 
between the results from observational data and the predic- 
tions from semi-analytic models of galaxy formation imple- 
mented within a high resolution simulation of the evolution 
of dark matter in a 'concordance' ACDM cosmology. 

In Section [2] we present our observational sample, ex- 
tracted from the latest Sloan Digital Sky Survey (SDSS) 
spectroscopic data. In Section [3l we explain how we obtain 
estimates of SFR and Z for our sample. In Section U we 
study the dependence of the M^-Z relation on SFR, and the 
dependence of the SFR-Z and sSFR-Z relations on M,. In 
Section [S] we discuss the latest versio n of the semi-an alytic 
galaxy formation code L- Galaxies (|Guo et al.ll20lj ). im- 
plement ed on the Millennium-II da rk matter N-body sim- 
ulation (jBovlan-Kolchin et al.ll2009l ). and describe how we 
select samples of simulated galaxies for comparison with the 
observational data. In Section [6] we present these compar- 
isons and show how trends in the relations between stellar 
mass, star formation rate and metallicity can be understood 
in terms of the prescriptions for gas accretion, star forma- 
tion, supernova and AGN feedback that are implemented 
in the model. In Section [71 we discuss the viability of dif- 
ferent physical mechanisms in regulating the metallicity of 
galajdes. Finally, in Section [SI we summarize our results. 



2 THE OBSERVATIONAL SAMPLE 

The sample of galaxies analysed in this paper is drawn from 
the Sloan Digital Sky Survey MPA-JHU Data Release 7 
catalogue (hereafter, SDSS-DR7]J. This catalogue contains 
~ 900, 000 galaxies with available spectra. The sample cuts 
used here are the same as those adopted bv lTremonti et al.l 
(|2004) . who investigated the M,-Z relation using galaxies 
from the SDSS Data Release 2. 

First, we remove all duplicate spectra from the cata- 
logue, reducing it by ~ 3.2 per cent and leaving 898,302 
galaxies. Then, following iTremonti et al.l (2004), we take 
only galaxies with reliable spectroscopic redshifts within the 
range 0.005 < z < 0.25. We then remove all gala:xies whose 
fibre-to-total light ratio is less than 0.1. This is defined as the 
ratio of the flux given by the SDSS fibre magnitude to that 
given by the SDSS model magnitude, in the r-band. This cut 
eliminates galaxies with metallicity measurements that are 
heavily biased towards the inner regions. We have checked 
that increasing the minimum fi bre covering fraction to 0.35 
(following recommendations bv lKewlev fc Ellisonll2008l ). or 
raising the minimum redshift, does not affect any of the main 
results presented in this paper. 

We also make cuts to the signal-to-noise ratio (SNR) of 
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Sample cut 



No. of objects removed 



Removing duplicates 


29,250 


Redshift cut 


137,087 


Fibre-to-total light cut 


34,203 


SNR cut 


372,438 


o-(mz, H<5a,D„4000) cut 


16,027 


AGN cut 


98,635 


M*, Z confidence cut 


62,841 


Tl cleansing 


56,580 


T2 cleansing 


64,264 



Table 1. The number of objects removed by each sample cut, in 
the order they were applied. The final number of objects in Sam- 
ple Tl and T2 is 120,491 and 112,797, respectively. The difference 
in final size between the two samples is due to the different re- 
quirements for estimating SFR and Z (see text). 93,971 galaxies 
are common to both samples. 



some of the key emission lines required to estimate metal- 
licity, ensuring that SNR(Ha, H^, [Nii]A6584) > 5. Again, 
we have checked that raising this threshol d to SNR ^ 10 
does n ot change our main results. Following iTremonti et al.l 
I1200J), we also make cuts on the accuracy of some addi- 
tional parameters that were used to estimate stellar masses 
in their original analysis. All galaxies for which airn^) < 0.15 
sinh"^(mag), a(RSA) < 2.5k, and cr(D„4000) < 0.1 are re- 
moved. This is done purely to achie ve consistency with th e 
original sample selection criteria of ITremonti et al.l (|2004I 1. 
The stellar masses that we use for our current DR7 analysis 
are derived using u, g, r, i, z SDSS photometry. 

Of those galaxies for which SNR([Oni]A5007) > 3, 
AG N hosts were rem o ved fol lowing the prescription given 
by iKaufTmann et al.l (l2003cll for defining AGN in the 
[Baldwin, Phillips fc TerlevichI l|l98ll ) (BPT) diagram: 



log([Oni]/H/3) > 0.61/log([Nn]/Ha) - 0.05 -f 1.3 . (1) 

For galaxies with SNR([Oin]A5007) < 3 , only those with 
log([Nn]A6584/HQ:) < —0.4 were retained, thus removing 
low-ionization AGN from the sample. 

Finally, a cut to the derived values of M, and Z was 
made, based on the 'confidence' with which they were esti- 
mated from fits to synthetic spectr a or HII region models 
using Cloudy ("Ferland et all 119981 ) (see Section [33]) . The 
\a spread in the likelihood distribution of the best-fitting 
model must be less than 0.2 dex in both quantities for the 
galaxy to remain in the sample. 

The number of objects removed by these cuts, including 
final cleansing (see Sections 13.11 and 13. 2p . is given in Table 
[l] Details of our two final samples, when binned by M« and 
SFR, are given in Tabled 



3 ESTIMATION OF METALLICITY AND 
STAR FORMATION RATE 

After creating our base sample of 177,071 emission line 
galaxies, we implement two different procedures for estimat- 
ing metallicity and star formation rate. 

The first procedure a i ms to replicate the methods out- 
lined in iMannucci et al.l (|2010l ). and uses the Ha-based 
method of iKennicutd ( 19981 ) for deriving SFR, and two of 



Sample Tl 


Minimum 


Median 


Maximum 


log(A/.) 

log(SFR) 

12 -h log(0/H) 

Redshift 


8.925 

-2.025 

8.72 

0.030 


10.125 

-0.525 

9.01 

0.077 


11.175 
0.975 
9.14 
0.179 


No. of galaxies 


120,491 






Sample T2 


Minimum 


Median 


Maximum 


log(A/.) 

log(SFR) 

12 -f log(0/H) 

Redshift 


8.625 

-1.125 

8.48 

0.019 


9.975 
0.075 
8.88 
0.068 


11.175 
1.575 
9.13 
0.174 


No. of galaxies 


112,797 







Table 2. Details of the parameter space coverage of our two 
observational samples. These values are given for data binned by 
Aft and SFR, considering only bins containing ^ 50 galaxies. 



the str ong line ratio calibrations described bv lMaiolino et al.l 
(12003) for deriving metallicity. We will refer to the list of 
metallicities, stellar masses and star formation rates gener- 
ated using this procedure as Sample Tl. 

The second procedure uses the values for SFR and Z 
provided by the SDSS-DR7 online catalogue. These esti- 
mates utilize emission line fluxes and Bayesian techniques 
implemented in a library of model galaxies that were created 
usin g a combination of stellar population synthesis models 



from 



Bruzual fc Charloti ll2003h and HII region models from 



ICharlot fc Longhettil ( 200ll ). The Bayesian estimation tech- 



nique provides a way of deriving robust errors for each of 
the derived physical parameters. We will refer to this list of 
estimates as Sample T2. Comparison of the two procedures 
provides insight into how sensitive the derived relations be- 
tween M», SFR and Z are to uncertainties in our SFR and 
metallicity estimates. 

We note that oxygen abundance is used as a proxy for 
global gas-phase metallicity throughout this work. We ex- 
press metallicity in terms of the number of oxygen atoms 
to hydrogen atoms in the gas component of a galaxy, nor- 
malised to the dimensionless quantity Z = 12 + log(0/H) . 
The current determination of t he solar oxygen abundance 
in these units is Z q = 8.69 (JAUende Prieto et al.l 1200 ll : 
lAsplund et~aLll2009l ). 



3.1 Sample Tl 

Following IMannucci et al.l (|2010|), total stellar masses are 
taken directly from the SDSS-DR7 catalogue, and corrected 
from a Kroupa to a Chabrier IMF by dividing by a factor 
of 1.06. Cleansing the sample of galaxies with: a) uncertain 
estimates of M, in the catalogue, b) emission lines that have 
too low signal-to-noise (SNR) for accurate Z estimates, c) 
estimates of Z from the two diagnostics used that differ by 
more than 0.25 dex (see below), reduces the sample to a 
total of 120,491 galaxies. 

SFRs were measured from the Ha emission line flux, 
corrected for dust extinction. The corrected Ha flux from 
each galaxy was converted to luminosity using L = Fc x 
A-kD^ X lO"^'^, where the distance D to the source was de- 
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rived from the galaxy's redshift, and 10~^^ is the factor used 
to normahse SDSS fluxes to units of ergs s~^ cm"^. Star for- 
mat ion_£ate_was then determined using the fixed conversion 
from lKennicuti j 19981 ). 



SFR = 7.9 X 10"*^ L(Hq) [ergs s" 



(2) 



which assumes a case B recombination, an electron temper- 
ature of 10* K and a Salpeter IMF. This value was then also 
corrected to a Chabrier IMF by dividing by 1.7. 

The observed Ha flux is corrected for external dust ex- 
tinction using 



iKewlev &: Dopital (|2002l ). The resulting co mbined calibra- 
tions are given by eqn. 1 and table 4 in the lMaiolino et al.l 



paper. 



For our Sample Tl, we follow iMannucci et al.l (|201Cll ) 
by taking the average of the metallicities given by the 
[Nll]A6584/Ha calibration and the R23 calibration as the fi- 
nal metallicity estimate for each galajcylfl We a l so cor rected 
all line fluxes for dust, following ICardelli et al.l ([l989). This 
lowers the metallicities estimated via R23 by ~ 0.03 dex at 
the highest masses, but makes very little difference to those 
estimated via [Nii]A6584/Hq: because the two lines involved 
are of very similar wavelengths. 



Fc(Ha)=Fob.(HQ)e+"H- 



(3) 



where the Ha optical depth rua can be determin ed using 
the w avelength-independent relation A\ — I.OSGta (|Calzettil 
[l99J). An estimate of Aua can be obtained from the Balmer 
decrement B = Fobs (Ha) /Fobs (H/3) and the following equa- 
tion: 



Ahc = -2.5 log 



B 



kuo 



kua — fci; 



(4) 



where 2.86 is the intrinsic Balmer decrement for a case B re- 
combination with electron de nsity Ue = 100cm ~^ and elec- 
tron temperature T^ = lO^K (|Osterbrocklll989l ). These val- 
ues are roughly appropriate for local sta, r forming galax- 
ies (JKcwlcv ct al.' 2001: 'iz otov et aJ] I2OO6I : ILIu et aLll20oi : 
IPilvu gin & Mattsson 2011) . The k- valu es used are /chq = 
2.468 and fcH/3 = 3.631 (|Calzetti|[200il 'l. and represent the 
ratio of the extinction at a given wavelength to the in- 
trinsic colour excess in the b-band r elative to the v-band 
(kx = Ax/E{B - V)^) (|Calzettilll993 ). 

It should be noted that these SFRs pertain to the re- 
gion of the galaxy falling within the 3" diameter SDSS fibre 
aperture, which probes the inner ^1-9 kpc of the galaxies in 
our samples. This method was used for Sample Tl in order 
to follow the procedure adopted by IMannucci et al.l (|20ld ') 
as closely as possible. The method used in Sample T2 in- 
stead estimates the total-SFR using the SDSS photometry 
to correct for the missing star formation in the outer regions 
of the galaxy (see Section [3. 2|) . 

Metallicity was calculated for Sample Tl using two 
of the strong line diagnostics calibrated by iMaiolino et al.l 
l|2008r . Such diagnostics are often used when direct measure- 
ments of the electron temperature Te in the Hll regions of a 
galaxy are not possible, or at high metallicities {Z > 8.35), 
where they no longer provide an accurate estimate of Z due 
to temperature fluctuations within individual Hll regions 
and a cross the whole galaxy (|Stasiriskall2005l : IMaiolino et al.l 
120081 ) (this is the case for essentially all the galaxies in our 
sample). In such cases, either an empirical method, utilis- 
ing other galaxies with measured Te metallicities, or a the- 
oretical method, utili sing purely theoretical p hotoionisation 
mo dels, can be used (JKewlev fc EllisonI [2003) . 

IMaiolino et al.l ((20081) derived their strong line diagnos- 
tics using a combination of empirical and theoretical meth- 
ods. 259 local gal axies of Z < 8 . 35, wi th Ts derived metallic- 
ities compiled by iNagao et al.l (|200y) were used, combined 
with 22,482 SDSS-DR4 galaxies of Z > 8.4, with metal- 
licities derived using the photoionisation model outlined in 



3.2 Sample T2 

Our second sample. Sample T2, utilises the stellar masses, 
SFRs and metallicities given in the online SDSS MPA-JHU 
DR7 catalogue. Values have been nominally corrected from 
a Kroupa to a Chabrier IMF to match Sample Tl. Total 
stellar masses and star formation rates are used, whereas 
metallicities are calculated using emission line fluxes that 
fall within the spectral fibre. Cleansing the sample of galax- 
ies without robust estimates for M, , SFR or Z from the 
catalogue reduced the sample to 112,797 galaxies. We note 
that the stellar masses in the SDSS-DR7 catalogue are based 
o n fits to photomet r ic data rather than to Lick indices, as 
in lKauffmann et al.l (|2003al '). 

Total-SFRs are ca lc ulated using the technique described 
by iBrinchmann et al.l (|200J), wit h im provements to the 
aperture corrections as detailed bv lSalim et al, (2007). This 
method for estimating star formation rates is based on 
fitting to a grid of photoionisation mo dels derived from 
the Cloudy code jFerland et al.l 1 19981 ') . as detailed by 
ICharlot fc Longhettil ( 200ll ). It thus accounts for the fact 
that the Ha-to-SFR conversion factor will depend on metal- 
licity, and it also allows the dust-free value of the Balmer 
decrement to differ from the 'standard' case B value (see 
IBrinchmann et al.l |2004| for a more extensive discussion) . 
The use of total- rather than fibre-SFRs is the main reason 
why the median star formation rate of galaxies in Sample T2 
is 0.6 dex higher than it is in Sample Tl (see Table 2). The 
most significant shift is seen at low masses, as these galaxies 
tend to have low redshifts and larger apparent sizes, and are 
therefore more extended with respect to the fibre aperture 
than more distant galaxies. 

Metallicities are calculated using the same grid of pho- 
toionisation models, by finding the model that best matches 
the observed fiuxes of the most prominent optical emission 
lines ([On], Rfi, [Qui], Ha, [Nil] and [Sii]). We refer the 
reader to Appendix A for a brief discussion on the merits 
of using this Bayesian technique over simpler emission line 
ratios when studying local, high-Z, star-forming samples. 



^ A modified version of the R23 ratio was used for Sample Tl 
wliicli does not require the [Olll]A4959 Une: R23 = ([Oll]A3727 + 
3.1 ■ [Oin]A5007)/H/3. This version provides a more robust diag- 
nostic for metaUicity when the SNR of the oxygen lines is low (F. 
Mannucci, private communication). 
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Figure 1. The M*-Z relation for Sample T2. The galaxy distri- 
bution is shown in grey, and a fit to this data is given (black solid 
line). The mean metallicity in bins of 0.15 dex in stellar mass 
is plotted as black diamonds. The Icr dispersion in Z about the 
mean is shown as dotted lines. Also shown is the lTremonti et al.l 
||2004| ') fit to data drawn from the SDSS-DR2 (dashed red line). 
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Figure 2. The distribution of residuals about Z for the M,- 
Z relation shown in Fig. [T] The mean dispersion of ~ 0.10 dex 
is the same as that reported by iTremonti et al.l 1 120041 ) for their 
SDSS-DR2 sample. 



4 OBSERVATIONAL RESULTS 

In this section, we examine whether the M,-Z relation ex- 
hibits additional dependences on SFR. The region of the pa- 
rameter space of interest for our analysis is 8.6 < log(M, ) < 
11.2, -2.0 < log(SFR) < 1.6 and 8.5 < Z < 9.2, as this 
covers 98 per cent of the galaxies in both our observational 
samples. 



4.1 The M*-Z relation 

Firstly, we introduce the basic M*-Z relation for Sample T2, 
before analysing any SFR-dependence. In t his form, the re- 
lation can be seen as an 'update' to that of ITremonti et al.l 
I1200J), who analysed a sample containing half as many 
galaxies. The M^-Z relation for Sample T2 is shown in Fig. 
[l] Individual galaxies are shown in grey, with a 3rd order 



polynomial fit to the whole population shown as a solid black 
line, and given by the following equation: 



26.6864 - 6.63995a: + 0.768653a:^ - 0.0282147a;-' 



(5) 



where Z = 12 -^ log(0/H) and x = log(M*/M0). A fit 
to the lTremonti et al.l (|2004l ) mass-metallicity relation (red 
dashed line) is plotted for comparison. The standard devia- 
tion about the best fit from residuals is 0.102 dex, as shown 
by Fig. [2] This is the same as the dispersion calculated for 
the original ITremonti et al.l (|2004| ) sample. 

Our new fit indicates a somewhat more lin- 
ear relation between mass and metallicity below ~ 
lO^'^M©. Above this mass, the relation fiattens and 
there is a hint of a turnover at the very high- 
est masses. Turnovers in the M«-Z rel ation at high 
mass are not uncommon in t he literature (jZaritskv et al.l 
19941: iKewlev fc Dopital l2002l: iKobuhiickv fc KewlevI l2004l : 



McGaughlll99ll : IPettini fc Page]||2004l ). however they have 
not been widely discussed in terms of their physical signifi- 
cance or with reference to star formation rates. 



4.2 The M,-Z relation, as a function of SFR 

In order to study the dependence of the M,-Z rela t ion o n 
SFR, we follow a similar approach to lMannucci et al.l l|2010l ). 
binning galaxies by M, and SFR, and calculating the mean 
metallicity in each bin. Bins are of width 0.15 dex in both 
dimensions, and only those which contain Js 50 galaxies are 
plotted. Fig.EJshows this M,-Z relation for Sample Tl (left 
panel) and Sample T2 (right panel). The data are coloured 
by star formation rate, as are fits to the relation at a series 
of fixed SFRs (solid lines) . 

The left panel of Fig. [5] clearly shows the result previ- 
ously described by iMannucci et al.l (|2010l ): there is an in- 
crease in metallicity with increasing mass, but also a clear 
and ordered dependence of metallicity on SFR at fixed mass. 
Metallicity depends strongly on SFR at low masses, but is 
virtually independent of SFR at high masses in this sam- 
ple. The 1(7 spread about the median Z f or Sample Tl is 
0.07 d ex, compared to 0.08 dex reported bv IMannucci et al.l 
l|2010t ). 

The M^-Z relation for Sample T2 is somewhat differ- 
ent. It extends down to lower stellar masses and metallicities, 
and has a higher median star formation rate. Additionally, 
rather than an ordered segregation of the relation for fixed 
SFR, Sample T2 exhibits a 'twist', where low-SFR galaxies 
change from being the most metal-rich at low masses, to the 
least metal-rich at high masses. Sample T2, therefore, ex- 
hibits a SFR-dependence at both low and at high mass, with 
a broad transition region between the two regimes around 
M* ^ 10^'^'^Mq. We note that we obtain the same features 
for our observational sampl es when using the sam ple se- 
lection criteria described bv IMannucci et al.l (|2010l ). rather 
than those outlined in Section [21 This indicates that the 
choice of sample cuts does not significantly affect our find- 
ings. 

The extension to lower stellar masses seen for Sample 
T2 is mainly due to the removal of galaxies from Sample Tl 
during cleansing. Low-redshift galaxies [z < 0.03) are ex- 
cluded from Sample Tl because the [Oll]A3727 line (required 
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Figure 3. The M^-Z relation for Sample Tl (left panel) and Sample T2 (right panel). Filled circles show the mean metallicities of 
galaxies binned by Af, and SFR, for bins containing ^ 50 galaxies. Binned data is coloured by SFR, as are the fits at fixed SFRs (solid 
lines), plotted for log(SFR) = -1.875, -1.575, -1.275, -0.975, -0.675, -0.375, -0.075, 0.225, 0.525 for Sample Tl, and log(SFR) = -0.975, 
-0.675, -0.375, -0.075, 0.225, 0.525, 0.825 for Sample T2. The dependence of Z on SFR for the two samples is clearly different. For 
example. Sample T2 exhibits a clear dependence at high-mass. 



for estimating metallicity using the R23 ratio) is not mea- 
surable. These galaxies can remain in Sample T2, extending 
the M,-Z relation down a further ~ 0.3 dex in stellar mass. 

The shift in SFR seen in Sample T2 is explained by a 
combination of factors. First, total- rather than fibre-SFRs 
are used (see Section 13. 2|) . Second, the simpler Ha-based 
method used for Sample Tl yields lower SFR estimates for 
high-M« galaxies of given Ha luminosity. This is because 
the conversion factor in E qn. [2] is dependent on m etallic- 
ity, and hence stellar mass. lBrinchmann et al.l ([2004!) report 
that this Kennicutt value is most accurate for galaxies with 
stellar masses ~ lO^'^M©, and will underestimate the SFR 
in more massive galaxies. Note, however, that the intrin- 
sic Balmer decrement for case B recombination - the case 
B ratio - is also metallicity- and mass-dependent. A fixed 
case B ratio can over-estimate the attenuation due to dust 
by up to ^ 0.5 mag for the most massive, metal-rich galax- 
ies iJBrinchniann et al.ll2004l ) . This will counteract somewhat 
the underestimate in SFR from a fixed conversion factor for 
higher mass galaxies. None the less, the use of fibre-SFRs 
and a fixed conversion factor combine to shift the median 
value of SFR down from log(SFR) = 0.075 Me/yr for Sam- 
ple T2 to log(SFR) = -0.525 Mo/yr for Sample Tl. 

The difference in the dependence of metallicity on SFR 
at high stellar masses in Sample T2 is mainly attributable 
to the metallicity derivation method chosen. Fig. [4] shows a 
comparison of the values of Z obtained from the two meth- 
ods used. Only the 93,971 galaxies that are present in both 
samples are included, and these are re-binned by M* and the 
Sample T2 value of SFR. We can see that for metallicities 
below Z ~ 9.1, the strong line ratio method used for Sam- 
ple Tl yields systematically higher values of Z compared to 
the Bayesian method used for Sample T2. This difference is 
also larger for galaxies with lower SFR. A similar effect is 
seen when binning galaxies by the Ha-based SFRs used for 
Sample Tl. 
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Figure 4. A comparison of the oxygen abundances obtained 
from the two metallicity diagnostics used in this work. Data is 
binned by M* and the estimation of SFR described in Section 
13.21 Points are co loured by SFR. The x- axis represents the tech- 
nique outlined b vlMannucci et al.l 1I2OICI) . using strong line ratio 
calibrations from lMaiolino et al.l ||2008|). The v-axis re presents the 
Bayesian technique outlined bv lTremonti et al.l ll2004f) . There is a 
clear and systematic over-estimation of Z from the former method 
relative to the latter. This discrepancy is also more significant for 
low-SFR galaxies. 



It is not straightforward to determine which of these two 
methods is most accurate at estimating oxygen abundance. 
The calculation of indirect gas-phase metallicities is fraught 
with complications, and a full discussion on the merits of dif- 
ferent methods is well beyond the scope of this work (how- 
ever, see Appendix A for a brief discussion). Noting that 
the properties of Sam ple Tl are already well described by 
iMannucci et al.l (|201GI ) using their data, we choose to drop 
further analysis of Sample Tl and focus on Sample T2 in 
the rest of this work. 
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Figure 5. The SFR-Z relation (left panel), {SFR/M,)-Z relation (middle panel) and D„4000-Z relation (right panel) for Sample T2. 
Fits to the data for four fixed masses are shown (solid lines). The drop in Z with increasing star formation rate (decreasing Dn4000) 
is seen in all three relations. A drop in Z with decreasing star formation rate (increasing D,i4000) is also seen for high-mass galaxies, 
representing the same effect seen at high masses in Fig. [S] 



4.3 The SFR-Z and sSFR-Z relations, as a 
function of AI, 

The SFK-Z relation for Sample T2 is shown in the left panel 
of Fig.O Points are coloured by stellar mass, and fits to the 
relation at four fixed stellar masses are also shown (solid 
lines) . 

We can see the drop in Z with increasing SFR for low- 
mass galaxies, however Sample T2 also exhibits a downturn 
in metallicity with decreasing SFR at high stellar masses. 
These plots show the same SFR-dependences seen in the 
Mt-Z relation, but from a nother projection onto the M«- 
SFR-Z space. We note that lLara-Lopez et al.l ([20101) carried 
out an independent investigation of the relation between M«, 
SFR and Z using the SDSS-DR7 online catalogue. They also 
noticed a significant downturn in metallicity with decreasing 
SFR in their sample, although they did not study this effect 
as a function of stellar mass. 

Finally, we note that it is actually more sensible to study 
the dependence of the Mt-Z relation on the specific star 
formation rate (SFR/M*, hereafter sSFRe), rather than the 
star formation rate. A dwarf galaxy with a star formation 
rate of 1 M0yr~^ is a much more 'active' system than a 
giant elliptical galaxy forming stars at the same rate. The 
specific star formation rate, on the other hand, is a nor- 
malised quantity, and provides a measure of the present- 
to-past-averaged star formation rate of the galaxy. Another 
related quantity is the 4000- Angstrom break strength. This 
is characterised by the D„4000 index, the average fiux from 
two narrow bands on either side of the break (3850-3950A 
and 4000-4100A). sSFR estimates based on Ha fiux are a 
good measure of the instantaneous star formation rate in a 
galaxy, whereas D„4000 is more sensitive to stars that have 
formed over timescales of a few hundred million years to a 
gigayear. 

The sSFRe-^' and D„4000-2' relations for Sample T2 
are presented in the middle and right panels of Fig. [5] The 
data has been binned by M, and the sSFR indicator in ques- 
tion. The reader should note that small values of D„4000 
correspond to high values of sSFR. 

It is encouraging that both these relations exhibit sim- 
ilar trends. Metallicity decreases as a function of sSFR 
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Figure 6. Top panel: The mean dispersion (solid line) in Z 
for different projections onto the M*-SFR-Z space for Sample 
T2. The la scatter in the mean dispersion is also shown (dotted 
lines). Projections onto the plane are defined by a, given in Eqn. 
\E\ The projection of least scatter for Sample T2 is o = 0.19 (for 
Sample Tl, this is a = 0.26). Bottom panel: The projection 
corresponding to the optimum value of a. Only a slight improve- 
ment in the dispersion is obtained from this projection compared 
to the Mt-Z relation. Points are coloured by SFR, as in Fig. \3\ 



for low-mass galaxies and increases as a function of sSFR 
for high-mass galaxies. When noting that D„4000 (an ab- 
sorption feature) and sSFRe (computed from emission line 
fiuxes) are independent quantities, the fact that the same 
trends with metallicity are seen for both indicates that our 
results are likely to be robust. 
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4.4 Projection of least scatter 

If both M, and SFR are correlated with metallicity, then 
a linear combination of the two may provide a tighter re- 
lation with metallicity than the traditi o nal M ,-Z relation. 
This was explored by iMannucci et alj (|2010|), who calcu- 
lated the scatter in median metallicity around their FMR 
for a series of projections onto the M*-SFR-Z space, fixing 
Z as a principle axis. We modify this method slightly to find 
the mean dispersion in Z of our binned d ata for the same 
2d projections. Following IMannucci et al.l (2010), the linear 
combination of M, and SFR used is denoted by ^a , where 



Ma =log(Af.)-alog(SFR) 



(6) 



The free parameter a defines the projection, and can be 
varied to shift from the M^-Z relation (a — 0) to sSFR~^- 
Z {a = 1). The corresponding dispersion function for the 
binned data of Sample T2 can be seen in the top panel of 
Fig. [S] To obtain this function, the spread in Z was calcu- 
lated in 0.1 dex bins in fia, and the mean spread for each 
projection found. These mean dispersions from a = to 
1 were then fit by a third order polynomial to provide the 
solid line shown. The la spread in the dispersions calculated 
for each projection was also found and fit in the same way 
(dotted lines). 

It is clear that the M,-Z relation {a = 0) is not the 
optimum projection for Sample T2. However, the decrease 
in scatter obtained when using the optimum projection 
(q — 0.19) is only slight (~ 0.04 dex). This can be seen 
by comparing the spread in the M,-Z relation in the right 
panel of Fig. |3] with that in the bottom panel of Fig. |6] It 
is also interesting to note that when using fibre-SFRs, the 
projection of least scatter for Sample T2 drops to a = 0.03, 
very close to the M«-Z relation. This is because the spread in 
star formation rates at low-M* is reduced due to the under- 
estimation of SFR for nearby galaxies. 

The improvement obtained for our Sample Tl is more 
significant. A projection of a = 0.26 reduces the dispersion 
in Z by ~ 0.09 dex compared to the M^-Z relation. The rea- 
son for the lack of significant improvement seen for Sample 
T2 is the inverse dependence of Z on SFR seen at low and 
high masses. The 'u-shape' that the relation therefore forms 
in the 3d space makes it difficult to find a projection which 
reduces the overall scatter as much as in Sample Tl. 

We therefore conclude that, although some improve- 
ment to the scatter can be obtained by combining stellar 
mass and star formation rate in this way, the dispersion 
around the A'h-Z relation for Sample T2 is already rela- 
tively tight, and that the M,-Z relation is suitable for most 
purposes. 



5 MODEL SAMPLE 

In this section, we investigate the relationship between Af«, 
SFR and Z in a semi-analytic model of galaxy formation 
(SAM). This allows us to analyse the detailed evolution of 
individual galaxies as well as global relations, and to com- 
pare the models to observations. Analytic descriptions of 
physical processes can be self-consistently incorporated into 
semi- analytic models and then easily adapted, making them 
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Figure 7. The M^-Z relation for type and type 1 model galaxies 
at z = 0. Binned data (filled circles) are coloured by their SFR, 
as are the fixed-SFR fits at log(SFR) = -1.725, -1.425, -0.975, 
-0.675, -0.375, -0.075, 0.225, 0.525. The distinctive feature here 
is the high-mass turnover for low-SFR galaxies. This turnover 
causes the reversal of SFR-dependence from low- to high-mass in 
the model sample, which is also seen in our observational Sample 
T2. 



more fiexible than current SPH simulations. The models 
also provide large samples of galaxies with predicted stel- 
lar masses, metallicities and star formation rates, enabling 
detailed statistical comparisons with the observations to be 

made. 

The model used in our study is L-Galaxies (|Guo et al.l 

l2011f ). a SAM gr afted onto dark mat ter halo data from 
th e Millennium (ISprineel et al.l bOOSl ) and Millennium- 
II IjEovlan-Kolchin et al.ll2009f ) dark matter N-body simula- 
tions. The current model is able to reproduce the observed 
stellar mass and luminosity functions, as well as the TuUy- 
Fisher and mass-metallicity relations for present-day star- 
forming galaxies. It also includes updated analytical treat- 
ments for gas cooling and stripping, as well as supernova 
and AGN feedback. The processes important in regulating 
metallicity that are included in the model are SN-driven 
outflows with a fixed ejecta speed (this means that metals 
escape more readily from low-mass galaxies), and gas in- 
fall both in the form diffusely accreted pristine gas and gas 
that was enriched and ejected by the galaxy at an earlier 
stage. This reintegrated material is returned more quickly 
into galaxies residing in more massive DM haloes. A detailed 
description of the analytic re cipes used in L-Galaxies can 
be found in lGuo et al.l (|201ll ). The model assumes that the 
stellar IMF is the same everywhere, and at all epochs. 

Our z = model sample comprises 43,767 central galax- 
ies extracted from the Millennium Database provided by the 
German Astrophysical Virtual Observatory (GAVOJj. We 
use the catalogues generated by running the L-Galaxies 
code on the dark matter halo merger trees from the 
MiLLENNlUM-II simulation, which is able to resolve DM ha- 



available at; http://www.g-vo.org/Millennium 
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Figure 8. The SFR-Z relation (left panel), (SFK/Mt)-Z relation (middle panel) and {Mcoid/Mi,)-Z relation for the 2 = model sample. 
Fits to the data at four fixed stellar masses are shown (solid lines). The reversal in SFR-dependence from low- to high-masses evident in 
Fig.[7]is again seen. At fixed stellar mass, Z is seen to decrease with increasing SFR (sSFR) at low masses, but decrease with decreasing 
SFR (sSFR) at high masses. Similarly, there is also a clear anti-correlation between gas-to-stellar mass ratio and metallicity at fixed- Af* 
above log(M;.oi;j/M») -^ —1.0, but a direct correlation at lower gas-to-stellar mass ratios. These results are complimentary to those seen 
for our observational Sample T2 in Fig. \5\ 



los down to a halo resolution limit of 1.89 x 1O*M0. Galaxies 
were selected by stellar mass (8.6 ^ log(Af, ) ^ 11.2) and 
star formation rate (—2.0 ^ log(SFR) ^ 1.6) to span the 
same region of parameter sp ace as our observat ional sam- 
ples. A WMAPl cosmology IjSpereel et al.l l2003l ) with the 
following parameters is assumed: flm = 0.25, flh = 0.045, 
r^A = 0.75, Us = 1, (78 = 0.9 and Ho = 73 km s"^ Mpc"\ 

Both type and type 1 central galaxies are included 
in the model sample. Type galaxies are those lying at 
the centres of their dark matter haloes. These galaxies ac- 
crete material that cools from the surrounding halo, as well 
as satellite systems that sink to the centre of the halo as 
a result of dynamical friction. Type 1 galaxies are satellite 
systems embedded within a larger halo. In the current imple- 
mentation of L-Galaxies, these galaxies still accrete cold 
gas from their surrounding self-bound 'subhalo'. Over time, 
the dark matter and gas in the subhalo is tidally stripped 
and accretion rates decline. The Millennium Database also 
contains so-called type 2 satellite galaxies, which have lost 
their surrounding host DM haloes through tidal stripping. 
These galaxies generally have no ongoing star formation and 
are therefore not included in our analysis. 



? -1.0 
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Figure 9. The sSFR-(Mcoid/A'/*) relation for our z = model 
sample. The relation for type galaxies (blue line) and type 1 
galaxies (green line) is shown, as well as the mean (solid black line) 
and median (dashed black line) relation for the full sample (shown 
in grey). The ler spread around the mean is shown as dotted lines. 
The low-sSFR galaxies that form the high-mass turnover in the 
model Mf-Z relation are indicated by red circles. These galaxies 
tend to have lower-than-average gas-to-stellar mass ratios for their 
sSFR. 



6 MODEL RESULTS 

6.1 The relation between A/., SFR and Z at z = Q 

We bin the galaxies in our sample by Af* and SFR in the 
same way as was done for our observations. Only bins con- 
taining ^ 25 galaxies are included in the model plots - half 
the number required in the observational samples. This is 
done in order to expand the dynamic range in the plots, but 
none of our results depends on this choice. 

The total cold gas-phase metallicity for each bin is cal- 
culated. This is given by the ratio of mass in metals to total 



cold gas mass, normalised to the solar metal abundancqj: 
^coid = 9.0 + log(Afcoid,z/A^coid/0.02). 

The M,-Z relation for our 2 = model sample is shown 



•* A change in the value of Zq down to that preferred by 
lAsplund et al] ||2009| ') (Zq = 8.69) does not effect the slope of 
the relation or the distribution as a function of SFR (see also, 
iBertone. De Lucia fc Thomasll2007l ) . We therefore cho ose to stick 
with the value of Zq 
tency. 



: 9.0 used bv lGuo et al.l l|201ll ') for consis- 
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Figure 10. The evolution from z = 2 to in mass (top panels), 
metallicity (middle panels) and SFR (bottom panels) for three 
typical galaxies from the high-mass, high-Z sub-sample. Galaxies 
in this sub-sample tend to undergo gradual enrichment of their 
cold gas phase over time. The galaxy IDs for these galaxies from 
the Millennium Database are provided at the top of each panel. 



Figure 11. The evolution from 2; = 2 to in mass (top panels), 
metallicity (middle panels) and SFR (bottom panels) for three 
typical galaxies from the high-mass, low-Z sub-sample. These 
galaxies show the gradual dilution of the cold gas phase that is 
characteristic of the low-Z sub-sample. The Galaxy IDs for these 
galaxies from the Millennium Database are provided at the top 
of each panel. 



in Fig. [71 We see a positive correlation between stellar mass 
and metallicity, as well as a segregation of this relation by 
SFR. Interestingly, this segregation is similar to that seen 
in Sample T2, with a reversal in the SFR-dependence from 
low to high masses. Another way of stating this result is 
that there is a 'turnover' at high stellar masses in the M,-Z 
relation for low-SFR galaxies. Such a feature is hinted at in 
Sample T2, but we note that the number of real massive 
galaxies with emission lines that are strong enough to mea- 
sure metallicity and that are not dominated by AGN emis- 
sion is quite small. This may explain the relative weakness of 
the feature in the observations compared to the model. The 
stellar mass at which the turnover occurs is ~ lO^^'^M©, 
again in quite good agreement with what is seen in Sample 
T2. It should be noted that the same turnover is seen when 
only type galaxies are included in the model sample, and 
so is not a consequence of environmental effects. 

The SFK-Z and the sSFR-Z relations are shown in the 
left and middle panels of Fig. [S] The D„4000 index is not 
available for our model galaxies, so only sSFRe = SFR/M» 
is used. As in our observations, Z depends on star formation 
at fixed stellar mass. As in Sample T2, the trend for high- 
mass galaxies is the opposite to that seen at low masses. The 
metallicities of galaxies of mass > IO^^'^Mq decrease with 
decreasing SFR (sSFR). 

In the right panel of Fig. [8] we plot a new relation, that 
between the gas-to-stellar mass ratio (Mcoid/M,) and gas- 
phase metallicity for the same model galaxies. Because cold 
gas is the fuel for ongoing star formation in a galaxy, the gas- 
to-stellar mass ratio should correlate with the enrichment of 
the interstellar medium. Indeed, this is what is seen in our 
model sample. The average Mcoid/Af* decreases with stellar 
mass and with Zcoid for galaxies with stellar masses less 



than ~ 10^ Mq. At higher stellar masses, we again see 
a turnover towards low metallicities at low Mcom/M*. We 
conclude that the low-sSFR galaxies that contribute to the 
high-mass turnover in the M,-Z relation also tend to have 
lower Mcoid/M, than other galaxies of a similar mass. This 
is also seen in the sSFR-(Mcoid/Af*) relation, shown in Fig. 
[9l where these galaxies are indicated by red circles. 

Currently, there is only limited data available on the 
gas fractions of high-mass galaxies with low star forma- 
tion rates. However, with the ongoing development of 
the G ass (GALEX Arecibo SDSS Survey, ICatinella et all 
20101) and ColdGas s (CO Legacy Database for GASS, 



Saintonge et al.ll201ll ) programmes, we will soon have a sig- 



nificant number of gas-to-stellar mass ratios for high-mass 
SDSS galaxies with which to compare. It will be interest- 
ing to see if the relationship between Mcoid/M^, SFR and Z 
presented in Figs. [8] and [9] is also found in these observations. 



6.2 Metallicity evolution in model galaxies 

In this section, we study the origin of the turnover in the 
model M-t-Z relation seen in Fig. [T] We do this by splitting 
the high-mass end of the sample into low metallicity and 
high metallicity sub-populations and studying differences 
in their evolutionary histories. We extracted two high-mass 
(Af, > 10^°*Mq) sub-samples: a high-Z {Z > 9.2) sub- 
sample containing 134 galaxies, and a low-Z (Z ^ 9.0) sub- 
sample containing 136 galaxies. The mass, metallicity and 
SFR evolution of these two sub-samples was then compared. 
The L- Galaxies model tracks six distinct mass com- 
ponents of galaxies: stellar mass (in the form of a bulge and 
disc), black hole mass, cold gas mass (ISM), hot gas mass 
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(ICM), ejected gas mass (IGM) and halo stars (producing 
the intra-cluster hght). Mass and metals can pass between 
these components along pre-defined routes, depending on 
the processes taking place. The top three panels of Fig. 
1101 show the time evolution of these mass components for 
three representative galaxies from the model with high stel- 
lar masses and high metallicities. In the middle panels, we 
show the time evolution of the stellar and gas-phase metal- 
licities of the same galaxies. The bottom panels show the 
time evolution of their star formation rates. 

The left panels display a type galaxy in which the 
stellar mass (solid orange line) has been steadily increasing 
since redshift two. In this galaxy, the mass of cold gas (solid 
blue line) is always higher than the critical value required for 
star formation, Afcrit (dashed blue line)|j There is a steady, 
gradual increase in the metallicity of the stellar, cold gas 
and hot gas components. This is because stars are formed 
continuously, synthesising and distributing metals through- 
out the galaxy at a higher rate than the dilution due to 
the accretion of metal-poor gas. Around 64 per cent of the 
galaxies in our high-Z sub-sample have formation histories 
similar to this. 

The middle panels show a galaxy that first evolves in a 
similar way to the galaxy shown in the left panels. It is then 
accreted onto a more massive DM halo at 2 '^ 1.0, becoming 
a type 1 object, at which point gas and dark matter begin 
to be tidally stripped. Some of these type 1 galaxies, like the 
example shown in the middle panel, then exhibit a sharp in- 
crease in gas-phase metallicity. This is because gas accretion 
onto the galaxy is reduced, but star formation continues, and 
as a result, metals continue to be dispersed into ever decreas- 
ing volumes of hot and cold gas. In other type 1 galaxies in 
the high-Z sub-sample, the cold gas mass drops below Merit 
after being accreted, causing the cold gas metallicity to re- 
main constant thereafter, due to a shut-down in both star 
formation and galactic gas accretion. Type 1 galaxies with 
star formation histories of these forms constitute ~ 30 per 
cent of the high-^ sub-sample. 

The right panels of Fig.[TO]show a class of central galaxy 
that formed a bulge and a supermassive black hole at red- 
shifts greater than two. The black hole (solid black line) then 
grew steadily through so-called 'radio mode' accretion of hot 
gas from the surrounding halo. Such galaxies are not repre- 
sentative of the high-Z sub-sample as a whole (they comprise 
only 6 per cent of the sub-sample) and are almost certainly 
classified as AGN and so missing from our observational 
samples. Nevertheless, we think they are rather interesting. 
The galaxy makes it into our high-Z model sub-sample due 
to the fact that it accretes many, many satellites over time 
(4,141 since redshift two, compared to only 21 and 71 for 
the galaxies in the left and middle panels, respectively). The 
satellites bring in fresh gas, leading to the very 'bursty' star 
formation history seen in the bottom-right panel of Fig. 1101 
We note that this galaxy resides in the center of the 9th 
most massive DM halo in the whole Millennium-II simula- 
tion at z = 0, and so such a high rate of satellite accretion 
is perhaps unsurprising. 



® Note that the Too mre Q disc instability criterion for star for- 
mation l lToomrelll964 ) is ex pressed in L-Galaxies in terms of a 
critical cold gas mass llGuo et al.ll201ll . Section 3.4) 
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Figure 12. The M^^^^ig^-M^n relation for our 2 = model sam- 
ple. The full sample is shown in grey, with the mean relation 
plotted as a solid black line. The Icr spread around the mean is 
shown as dotted lines. Galaxies in the low-Z sub-sample are high- 
lighted by red circles. These galaxies clearly have higher black 
hole masses that expected from an extrapolation of the model 
relation, but lie nicel y along the observational relation derived by 
iHaring fc Ri^ ||2004|) . 



In conclusion, for the majority of galaxies in the high-Z 
sub-sample, the dominant process driving metallicity evo- 
lution is clearly a gradual enrichment of the gas phase due 
to continuous star formation. At these high masses, out- 
flows are inefficient at removing cold gas and metals from 
the galaxy. At these high SFRs, galactic infall rates are too 
low to dilute the ISM. Consequently, these galaxies become 
increasingly metal-rich with time. 

Fig. [TT] paints a rather different picture for the evolu- 
tion of galaxies in the low-Z sub-sample. All three galaxies 
have undergone dramatic drops in their cold gas masses, co- 
inciding with a merger event at some stage during the past 
ten gigayears. Mergers with mass ratio less than 10:1 are 
marked by black dots on the Mcom evolution line. Note that 
no such mergers occur for the three high-Z examples in Fig. 
1101 - only 7.5 per cent of galaxies in the high-Z sub-sample 
have undergone a significant merger over the last ten gi- 
gayears. A merger with mass ratio 3:1 or less is considered 
a major merger in the model, causing the destruction of 
the stellar and gas discs and the transfer of this material 
to the bulge of the descendant. We see from Fig. 1111 that 
not only major mergers cause the sudden drop in Afcoid. 
Gas-rich minor mergers are also effective at inducing star- 
bursts and the rapid growth of the central SMBH through 
'quasar mode' accretion. During such events, a black hole 
can grow by swallowing both cold gas and the smaller black 
hole of its companion. 92 per cent of the galaxies in the low-Z 
sub-sample have present-day black holes with masses greater 
than 10®"Mq, that were formed through this process. The 
remaining 8 per cent have either grown their black holes 
gradually through radio-mode accretion, or do not contain 
a central SMBH. 

We note that although large black holes are a feature of 
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Figure 13. The (SFR/A/t)-Z relation for our model sample (left panel) and observational Sample T2 (right panel) for galaxies of stellar 
mass > 10^^' ^Mq. Galaxies had been binned by sSFRe and Z^^i^, and bins are coloured by the mean black hole mass of the constituent 
galaxies. The number of galaxies in each bin is also given by white text. The low sSFR, low-Z population clearly have larger central 
SMBHs in the model sample. A similar dichotomy can be seen in Sample T2, though the distinction is not as clearly defined. The red 
boxes in the right panel mark the low-Z and high-Z regions shown in Fig. 1141 



almost all the galaxies in the low-Z sub-sample, they do not 
cause the low metallicities seen in these galaxies at z = 0. 
These are instead caused by a cessation in star formation due 
to the sudden drop in cold gas mass below Merit, followed 
by accretion of metal-poor gas. This galactic accretion is 
limited to a low rate by the suppression of cooling from radio 
mode AGN feedback, allowing it to continue for an extended 
amount of time without re- igniting star formation. The three 
galaxies in Fig. [11] therefore show that gradual dilution of 
the gas phase due to metal-poor infall of gas in the absence 
of star formation is the main process producing the low-Z 
sub-sample. 

This is an effect that is not seen in hydrodynamic sim- 
ulations of galaxy evo lutio n such as those carried out by 
iFinlator fc Pavel ([2008) and lDave. Qppenheimer k, Finlatorl 
1 2OIII ). In those models, galaxies quickly fall back into an 
equilibrium between their infall, outflow and star forma- 
tion rates after a perturbative event, whereby Minfaii = 
A/* + i\foutflow Instead, the inclusion of AGN feedback in 
our model enables galaxies to slowly accrete metal-poor gas 
for a number of gigayears without forming stars. 

We note that it remains to be seen whether the fraction 
of massive galaxies in the models with very little ongoing 
star formation and low metallicity reservoirs of cold gas is 
matched in observations. The low gas fractions of massive 
galaxies means that the gas is difficult to detect. Obtain- 
ing Hi maps of a large sample of such objects remains a 
significant observational challenge. 



6.3 Bulge and black hole masses 

As discussed in the previous section, one clear distinguish- 
ing feature of massive, low-Z galaxies in the model, aside 
from a low gas-to-stellar mass ratio, is the presence of a 
massive bulge and massive central black hole formed dur- 
ing a merging event. This is illustrated by the Afbuigc-MBu 
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Figure 14. The distribution of black hole masses in Sample T2 
for low-Z (blue histogram) and high-Z (red histogram) galaxies. 
These are the galaxies marked by the red boxes in Fig. 1131 The 
mean A/bh of each region is also show by vertical dashed lines. 
Low-Z, low-sSFR galaxies have more massive central black holes 
in Sample T2, as seen in the model sample. Shifting the bound- 
aries of the two regions somewhat does not affect this result. 



relation in Fig. 1121 where the galaxies in the low-Z sub- 
sample are highlighted by red circles. Interestingly, it is these 
gal axies that lie closest to the observational relation found 
by iHaring fc Rixl (|200J). A significant fraction of interme- 
diate to massive galaxies with smaller central black holes 
have likely had their bulges grown through secular processes 
(jShankar et al.ll201ll : Shankar, in prep.). Although massive 
central black holes are not the cause of the low metallicities 
in the low-Z sub-sample, they are an associated feature that 
we can look for in the currently available observational data. 
In the left panel of Fig. 1131 specific star formation rate 
is plotted against gas phase metallicity for all our high mass 
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model galaxies. Galaxies are binned by sSFRe and Z^oid, 
and the bins are coloured by the mean black hole mass of 
the galaxies in each bin. The number of galaxies in each bin 
is also indicated on the plot in white text. It is clear that 
low-sSFR, low-Z galaxies have the largest central black holes 
in the model. 

The right panel of Fig. [13] shows the same plot for an 
adapted version of our observational Sample T2 (see Ap- 
pendix B for details). Black hole masses were estimated via 
the measured stella r velocity dispersion u sing the Mbh-o" 
relation provided bv lTremaine et al.l l|2002h : 



log(MBH) = a + /?log(o-/(To) 



(7) 



where a = 8.13, ^ = 4.02, ctq = 200km/s and A/bh is in 
units of Mq. We can see that the observational data also 
contains a low-Z, high-MsH population, though the distinc- 
tion between this population and the majority of the sample 
is less clear than in the model. 

The dichotomy in the observations is more clearly seen 
in Fig. 1141 where the distribution of black hole masses 
is shown for galaxies contained within the two red boxes 
marked in Fig. 1131 The black hole population is clearly 
shifted to higher masses in the low-Z region (blue his- 
togram), compared to the high-Z region (red histogram). 
24 per cent of galaxies within the low-Z region have black 
holes of mass ^ lO^'^M©. This is only true for 2.7 per cent 
of galaxies in the high-Z region. A correlation between low 
metallicities, low-sSFRs and high black hole masses there- 



fore seems present in both our model and observational sam- 
ple. 

6.4 Evolution of the M,-Z relation out to z ~ 3 

In this section, we analyse the evolution of the M,-Z relation 
out to 2; ~ 3. Obse rvations have shown a clear evolution in 
the M^-Z relation (iMaiolino et aLlbOOSl) and M ,>-SFR rela- 
tion l|Noeske et al.ll2007al lbl: iKaiisawa et al.ll2010l ) with look- 
back time, and such evolution is required for the Mannucci 
FMR to remain fixed out to z ~ 2.5. 

In order to investigate this in L-Galaxies, three sup- 
plementary, identically selected samples of type and 1 
galaxies were extracted from the database at redshifts z — 
1.0, 2.0 and 3.1. These samples contain 63,745, 63,017 and 
50,558 galaxies, respectively, within the parameter space of 
interest. 

The evolution of the K'h-Z relation from z = 3.1 to 
in the model is plotted in Fig. 1151 The mean SFR at 
each redshift is also given in the top right of each panel. 
The SFR is seen to evolve strongly with redshift, drop- 
ping by ~ 0.9 dex from 2 = 3.1 to 0.0 at fixed m ass, in 
line with the drop ob served bv lNoeske et al.l (|2007al lbl) and 
IKaiisawa et al.l (|201G| ). However, there appears to be nearly 
no evolution in metallicity at fixed stellar mass at all in 
the model, contrary to observations. The present day M, - 
Z r elation agrees well with those of iTremonti et al.l (|200J) 
and lKewlev fc EUisonl (|2008l ) at z ~ 0, but the discrepancy 
with observations becomes increasingly pronounced towards 
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Figure 16. The evolution of four model galaxies from our 
2 = sample from redshift 2 to 0. Galaxies of final stellar mass 
log(M«) = 8.55, 9.49, 10.91 and 11.06 are shown in brown, or- 
ange, blue and red, respectively. The two higher mass examples 
are GId46078923 (blue) and GId76878142 (red) from the high- 
mass sub-samples shown in Figs.FlOland llll The full present day 
population is also shown in grey for comparison. We can see that 
galaxies tend to evolve along the 2 = relation at all masses, but 
that the nature of this evolution is different for different mass 
intervals. 



is seen observationally. iLee et al.l l|2006l ) have used a sample 
of 27 nearby dwarf galaxies to argue that the scatter in the 
observed local M,-Z relation remains roughly constant down 
to stellar masses of ~ lO^'^Mo. If confirmed with larger 
samples, this might suggest that, unlike metals, cold gas 
is not drive n very effectively out of l ow mass galaxies by 
supernovae (jMac Low &: Ferraralll999l ). 

The main point to take away from Fig.[T6]is that galax- 
ies evolve along the present-day M,-Z relation in the model, 
rather than from much lower metallicities, as suggested by 
observations. This suggests that it is indeed an overly rapid 
enrichment of the cold gas phase of galaxies before redshift 
three that is causing the lack of evolution thereafter in the 
model. 

One possible solution to the problem is that ac- 
creted gas is 'hung-up' in the atomic phase of the 
ISM for some time, before it is able to reach high 
enough densities to f orm the molecular clouds in 
which stars are formed IjGnedin. Tassis fc Kravtsovl |2009| : 



iGnedin fc Kravtsovllioiol ). Recently. iFuet all l|2010l ) imple- 
mented simplified prescriptions for the formation of molecu- 
lar gas in galaxies in the L-Galaxies code, to study scaling 
relations between gas and stars in local galaxies. The recipes 
have, however, not yet been fully applied to models running 
on the higher r esolution M i llenn ium-II simulation, or to 
the more recent I Guo et al.l (|201ll ) version of L-Galaxies, 
which contains stronger and more mass sensitive feedback. 
This will be the focus of future work. 



higher redshifts. This suggests that chemical enrichment in 
the model galaxies proceeded too rapidly at early times. 

In order to diagnose whether this hypothesis is correct, 
we plot the enrichment history from z — 2 to the present 
day of four representative galaxies from the z — model 
sample in Fig. 1161 The M*-Z relation for the entire 2 = 
sample is also plotted in grey. This is a reasonable test, as 
a large number of galaxies in the z — sample have main 
progenitors present in the higher redshift samples (the per- 
centage of galaxies that lie on the most massive progenitor 
(MMP) branch of a 2 = galaxy is indicated in the top left 
corner of each higher redshift panel in Fig. I15|l . 

We can again see in Fig. [16] the two types of behaviour 
at high stellar masses described in Section 16.21 with both 
■^coid and stellar mass gradually increasing with time for 
the blue (high-SFR) galaxy, and metallicity decreasing with 
time without any associated increase in stellar mass for the 
red (low-SFR) galaxy. At intermediate masses, stellar mass 
and metallicity increase together, punctuated by periods 
when Zcoid drops suddenly, due to an episode of enhanced 
gas accretion, such as a gas rich merger with a metal-poor 
satellite. 

At low masses, the metallicity evolution appears much 
more erratic. Drastic fluctuations in Z^om may occur when 
outflows become extreme enough to drive most of the gas out 
of the galaxy. This is the case for the lowest mass galaxies 
in the model, as a result of the SN feedback scheme that has 
been implemented in L-Galaxies. 

We can see from the grey points in Fig. [16] that the 
scatter in the M^-Z relation at z = increases towards 
lower masses as a consequence of this stochastic mode of 
evolution. We note that it is unclear whether this feature 



T DISCUSSION 

In Section |3] we have shown that our preferred observa- 
tional sample. Sample T2, exhibits a reversal in the depen- 
dence of metallicity on star formation rate from low to high 
stellar mass. At low masses, low-SFR galaxies have higher 
metallicities compared to more star forming galaxies. At 
high masses, they have lower metallicities. This observation 
alone could be explained by considering mass- dependent, 
metal-rich outflows. At lower masses (~ lO^M©), where out- 
flows are more effective, low-SFR galaxies produce relatively 
fewer SNe, disrupting the ISM less and blowing away met- 
als less efficiently (e.g. JMac Low fc Ferraralll99iJ ). At higher 
masses (> lO^'^'^Mo), where outflows become increasingly 
weak, low-SFR galaxies simply produce less metals and so 
under-enrich the ISM relative to more actively star forming 
galaxies. 

However, our model also shows a decrease in metallic- 
ity with stellar mass at flxed-SFR above ~ 10^"'* M©. Such 
a feature cannot be explained by mass-dependent outflows 
alone. If the turnover in the M^-Z relation is indeed real, 
then additional physical mechanisms must be at play. 

Our model points to metal-poor galactic infall at high- 
mass as an explanation. Those high-M* galaxies with low 
gas-phase metallicities are known to have undergone gradual 
dilution of their gas phases after a merger event which shut- 
down further star formation. The restriction in the amount 
of infall by AGN feedback allowed these galaxies to slowly 
dilute their ISM, without accreting enough gas for star for- 
mation to resume. Correlations between these galaxies and 
the high-M«, \ow-Z galaxies in our observational sample 
(namely, their large black hole masses) imply that such a 
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dilution process could also be involved in shaping the M,-Z 
relation in the real Universe (see also Appendix B). 

There are, however, two factors hampering this in- 
terpretation. First, the dependence of Z on SFR is itself 
strongly dependent on how these properties are measured 
(see Section 2)). Although a high- mass dependence is unde- 
niable in our Sample T2, it is not present in Sample Tl, at 
least not within the range of masses studied. 

Second, the recipes used to model physical processes in 
L-Galaxies, although up-to-date with current theory, are 
still rather crude, and this could be affecting the galaxy evo- 
lution seen in our model galaxies. For example, metals are 
assumed to fully mix with the ISM before galactic outflows 
are allowed to drive gas out of the galaxy. This means that 
it is the subsequent cessation of star formation in low-mass 
galaxies that is causing the relation between M, and Z in 
the model, rather than explicitly metal-rich outflows. Addi- 
tionally, and perhaps relatedly, there is a lack of evolution 
in the model Mt-Z relation, contrary to observations. 

Despite these two caveats, we believe our interpretation 
to be viable. On the observational side, the use of many 
emission line fluxes and SED fitting (as is done for many 
trusted stellar mass estimations) is likely to be a more ro- 
bust way of estimating metallicities than using individual 
emission line ratios (see Appendix A). On the modelling 
side, the problems outlined above are currently common 
to models in general. For example, the SPH sim ulation de- 
scribed bv lDave. Finlator fc Qppenheimeij (|201ll ') , which ex- 
hibits the same evolution of galaxies along the present-day 
Mt-Z relation, also has difficulties reproducing the observed 
metallicity evolution. Their favoured vzw model (which in- 
vokes momentum-driven winds and metal-rich outflows at 
all masses) shows an increase in Z of ~ 0.15 dex from z = 2 
to for galaxies of stellar mass ~ 1O^'^M0. Their fixed, low 
wind velocity sw model shows an evolution of ~ 0.1 dex 
(an amount similar to that seen in L-Galaxies). In com- 
parison, observations suggest an evol ution of ^ 0.5 dex a t 
the same mass from z = 2.3 to 0.07 (|Maiolino et al.l 120081 ). 
Further improvements to the observational determination of 
gas-phase metallicities, and the ongoing improvements in ac- 
curate modelling of galactic evolution, are necessary before 
significant progress can be made in this area. 



8 CONCLUSIONS 

We have shown that the gas-phase metallicities of galaxies 
are dependent on their star formation rate. This is also true 
at high masses, where highly star forming galaxies are seen 
to have higher metallicities - the opposite trend to that seen 
at low masses. However, despite this dependence, a projec- 
tion onto the A/*-SFR-Z space that combines M* and SFR 
does little to reduce the scatter in Z compared to the M,-Z 
relation. 

We also demonstrate the significance of metallicity 
derivation methods when assessing the relation between M, , 
SFR and Z. Strong line ratio diagnostics, such as R23, pro- 
vide significantly different metallicity estimates to Bayesian 
techniques which utilise emission line fluxes. These differ- 
ences are greater for low star forming galaxies. Although we 
believe that the Bayesian technique used for our Sample T2 
provides a more robust measurement of the global gas-phase 



metallicity in local, high-metallicity, star forming galaxies, 
it remains unclear to what extent this is so. 

In Section[Bl we show that a high-mass SFR-dependence 
is also present in our model sample. This is due to a turnover 
in the mass-metallicity relation, caused by a gradual dilution 
of the gas phase in some galaxies, triggered by a gas-rich 
merger which shuts down subsequent star formation without 
impeding further cooling. We have proposed that similarities 
between these low-sSFR model galaxies and those observed 
at z — 0, such as their larger-than-average black hole masses, 
leave open the possibility that such a process is also driving 
the SFR-dependence seen at high masses in real galaxies. 
If this is the case, then physical processes other than mass- 
dependent outflows must also be playing a part in regulating 
metallicity. Our model indicates metal-poor galactic infall as 
a likely candidate. 
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APPENDIX A 

As mentioned in Section 14.21 obtaining accurate estimates 
of the gas-phase metalhcity is not an easy process, as differ- 
ent diagnostics can provide very different results, ffowever, 
it should be emphasised that the difference seen in Fig. |4]is 
not simply due to the inherent discrepancies between Te, em- 
pirical and theoretical metallicity derivation methods. For 
the range of metallicities covered by our samples, both take 
their metallicities via theoretical methods. The difference is 
instead likely due to the use of either emission line ratios 
or emission line fluxes (with a Bayesian approach) in the 
analysis. 

On this point, it should be noted that theoretically de- 
rived strong line ratios are known to suffer from a number 
of problems, including d egeneracies (jKewlev fc Dopita 2002 : 
iKewlev fc Ellisonll2008l), sens i tivity to the ionisa tion param- 
eter (JKewlev fc Dopit'al 120021: [Erb et al.ll2006cl ') , saturation 
at high metallicities I Kewlev fc Dopital 20021: iLiang et al.l 



l2006l : lErb erZI l2006d : iKewlev fc EU ison '2008' ) and 

sistency with Te der i ved metallicities ( | Kcnnicutt et al.ll2003l : 
iBresolin et"ai]|2004l : ICarnett et "aIll2004l : IStasihskall2005l '). It 
is not clear to what extent the technique used for Sam- 
ple T2 suffers from these effects. There have been some 
concerns over the treatment of secondary nitrogen in the 
population synthesis mo dels used to produce the Sampl e 
T2 metaUicity estimates l|Liang et al.ll2006l : lYin et al.ll2007^ . 
Early indications show that accounting for this by exclud- 
ing the [Nil] (and [Sll]) lines from the Bayesian analysis 
does not change the high-mass SFR-dependence seen. Re- 
moving these lines only seems to have a significant effect at 
low-masses, strengthening the dependence of Z on SFR seen 
there0 

There are, however, two particular issues affecting 
the strong line ratio calibrations used for Sample Tf at 
high metallicity. Fir st, there is the binning of data by 
iMaiolino et al.l (|2008l 'l when calibrating the diagnostics used 
for our Sample Tf . A fit to unbinned data would have heav- 
ily biased their diagnostics against the lower metallicities 
crucial for high redshift studies such as theirs, due to the 
paucity of \ow-Z galaxies available in their present-day cal- 
ibration sample. However, this does mean that their fits are 
less precise at high metallicities, which is important for their 
application to local samples such as ours. Their R23 diagnos- 
tic over-predicts the average metallicity slightly compared 
to the average values obtained from the 'Kewlev fc Dopita 
2OO2I) model for Z > 9.0 (see fig. 5 in the Maioli no et al 



2008fl paper). This raises the metallicity estimated in this 



regime somewhat, despite the use of the [Nil] /Ha diagnostic 
to bring dow n the final va l ue. 

Second, iNagao et al.l l|2006l ) have already pointed out 
that their metallicity deri vation method - whi ch uses the 
same calibration sample as IMaiolino et al.l ()2008l l - may over 
estimate the gas met allicity at Z > 9. by a factor of ~ 0.1 
dex compared to the iTremonti et al.l (|2004l ) method. This 



is partly due to the bias towards selecting strongly [On] 
and [Olll] emitting galaxies in the sample causing the fit 
to be steeper at high metallicities. When considering the 
R23 diagnostic, the fact that H/? emission is also weaker 
in low star formation environments could explain why the 
discrepancy we see is more significant at low-SFR. 

Taking this into account, it is perhaps prudent to sug- 
gest that such ratios are not ideal for estimating metallici- 
ties for high-Z, local samples, where the availability of good 
spectroscopic data, including all optical emission lines, al- 
lows alternative techniques to be utilised. 



APPENDIX B 

The adaptation made to Sample T2 in Section [6.31 is to in- 
clude an additional set of 9,275 high-mass (> lO^^'^Mo) 
galaxies to the sample. These are the galaxies removed from 
the default Sample T2 due to having la uncertainties in 
their stellar mass estimates greater than 0.2 dex. The rea- 
son for this larger error seems to be due to large errors in the 
u-band magnitude measured, which is included in the esti- 
mation of stellar mass. This error has therefore likely prop- 
agated through to the confidence in the best fitting model 
during SED fitting, causing these galaxies to have a larger 
uncertainty in their M^, estimate. 

Fig. IB. II shows that, for masses greater than ~ 
IO^^'^Mq, these additional galaxies have only slightly 
greater la errors on their mass estimates than the default 
Sample T2 galaxies, and that their metallicity estimates are 
actually more than good enough for them to remain in the 
sample. We therefore choose to included these high-mass 
galaxies for our analysis of black hole mass. Their addition 
is not a necessary condition, but does provide more low- 
sSFR galaxies and makes the dichotomy seen in the right 
panel of Fig. 1131 clearer. 

It is interesting to note that when plotting the M^-Z 
relation for Sample T2 including these additional galaxies, 
a range of SFR-dependent turnovers can be clearly seen at 
high mass. This is show in Fig. IB. 21 When considering Fig. 
1131 Fig. [14] and Fig. IB. 21 together, it seems plausible that 
the low-Z, high-A/» galaxies with large SMBHs seen in our 
model sample are also present in the SDSS, and that pro- 
cesses other than outflows are at play in regulating their 
metallicities. 



° It should be noted that a larger number of double-peaked like- 
lihood distributions are produced when removing emission lines 
from the analysis. This can make it more difficult to determine 
the true metallicity. Further analysis of this new set of metallic- 
ity estimates is needed before more concrete statements can be 
made. 
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Figure B.l. The mean \a errors on the values of Af, (sohd hnes) 
and Z (dashed Unes) provided by the best fitting models for the 
default Sample T2 (thick lines), and for the additional galaxies 
described in Appendix B (thin lines). At high masses, the error 
in M« for the additional galaxies is only slightly greater than for 
the Sample T2 galaxies. The metallicity errors in the high-mass 
regime are well below the maximum acceptable error of ct = 0.2 
dex (indicated by the horizontal dotted line). 
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Figure B.2. The Mt-Z relation for Sample T2, when including 
the additional galaxies described in Appendix B. The shape of 
the overall relation is unchanged, but turnovers in the relation at 
fixed-SFR are now evident for wide a range of SFRs. This lends 
favour to the conclusion from our model sample, that processes 
other than outflows are regulating metallicity at high mass. 



